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ABSTRACT 

We present a schematic model for the formation of baryonic galactic halos and hot 
gas in the Local Group and the intergalactic medium. We follow the dynamics, chemical 
evolution, heat flow and gas flows of a hierarchy of scales, including: protogalactic clouds, 
galactic halos, and the Local Group itself. Within this hierarchy, the Galaxy is built via 
mergers of protogalactic fragments. Hot and cold gas components are distinguished, with 
star formation occurring in cold molecular cloud cores, while stellar winds, supernovae, 
and mergers convert cold gas into a hot intercloud medium. We find that early bursts 
of star formation lead to a large population of remnants (mostly white dwarfs), which 
would reside presently in the halo and contribute to the dark component observed in the 
microlensing experiments. The hot, metal- rich gas from early starbursts and merging 
evaporates from the clouds and is eventually incorporated into the intergalactic medium. 
The model thus suggests that most microlensing objects could be white dwarfs (m ~ 
O.SMq), which comprise a significant fraction of the halo mass. Furthermore, the Local 
Group could have a component of metal-rich hot gas similar to, although less than, 
that observed in larger clusters. We discuss the known constraints on such a scenario 
and show that all local observations can be satisfied with present data in this model. 
The most stringent constraint comes from the metallicity distribution in the halo. The 
best-fit model has a halo that is 40% baryonic, with an upper limit of 77%. Our model 
predicts that the hot intragroup gas has a total luminosity 1.5 x 10^° erg s~^, and a 
temperature of 0.26 keV, just at the margin of detectability. Improved X-ray data could 
provide a key constraint on any remnant component in the halo. 

Subject headings: dark matter — galaxies: evolution — galaxies: interactions — nuclear 
reactions, nucleosynthesis: abundances cosmology — dark matter, galactic evolution 



1. Introduction 

Recently there has been renewed interest in the 
nature of the dark matter in galaxy halos, motivated 
by the results of microlensing experiments. Observa- 
tions toward the Magellanic clouds (e.g., Alcock et al. 



1995 ; Aubourg et al. 1993) an d toward the Galactic 



bulge (e.g., Udalski et al. 1993 ) have detected gravita- 
tional microlensing and inferred the presence of dark, 
massive compact halo objects (MACHOs). Further- 
more, a recent binary detection in the direction of the 
LMC, along with average event durations of about 2.5 
months, may imply masses of order ~ 0.5Mq, sugges- 
tive of white dwarfs. 

At the same time. X-ray observations of clusters 
(Mushotzky 1992) and gr oups ( e .g., M ulchaey, Davis, 
Mushotzky, & Burstein |1996|, |1993|; Pildis, Breg- 



man, & Evrard [l995 ; Ponman, Bourner, Ebeling, & 



Bohringer 1996| ) have discovered a large amount of 
hot, metal-rich gas. Where it is observed, this gas is 
a substantial fraction of the baryonic mass-it is the 
dominant baryonic component of clusters, and is com- 
parable to the galactic component in groups. Indeed, 
it appears that most baryons in the universe are in 
the for of this hot X-ray emitting gas. 

The relatively high metallicity {Z ~ O.3Z0 for clus- 
ters, Z ~ 0.1^0 for groups) of this gas is impressive 
and demands that the material has undergone a sig- 
nificant amount of stellar processing during an ear- 
lier epoch of star formation. Such an epoch would 
also produce remnants, mostly white dwarfs. If this 
epoch is a general consequence of the formation of the 
bulge and halo of spirals like the Milky Way, as well 
as the ellipticals of rich clusters, it could account for 
the observed microlensing objects. 

Although white dwarfs are attractive MACHO 
candidates, there are important constraints o n suc h 
objects and their formation (Ryu, Olive, & Silk 1990 ). 
These include background light from the early evolu- 
tion, the present luminosity of the halo, and the metal 
and helium content of the disk and halo stars. While 
these place important constraints on model parame- 
ters, they do not rule out a significant white dwarf 
halo population, as we shall show. 

In modeling these galaxy aggregates and their hot 
gas components, one must account for the dependence 
of these systems on the morphology of the constituent 
galaxies. On the one hand, the hot gas in clusters and 
groups appears correlated with the luminosity of el- 
liptical and SO galaxies (Arnauld et al. 1992| ). This 



suggests that the hot gas arises from the violent merg- 
ing associated with these morphological types. On the 
other hand, it also seems well established (e.g.. Rich 
1990| ) that the morphology of the bulge and halo of 



spiral galaxies is quite similar to that of ellipticals 
and SO's. This suggests that the bulge and halo of 
spirals may have experienced a similar epoch of star 
formation and outflow during their formation. The 
difference in the morphologies may relate to the larger 
angular momentum or shallower gravitational poten- 
tial of spirals (Zurek et al. 1988 ), such that some gas 
survives halo formation, and settles afterwards into 
spiral arms. 

With this background in mind, we find that a likely, 
and perhaps inevitable, consequence of the forma- 
tion of the bulge and halo is the formation of a large 
remnant population in the Galactic halo, along with 
hot X-ray emitting gas in the Local Group and inter- 
galactic (i.e., extragroup) medium. While our model 
should be widely applicable, in this paper we concen- 
trate on the Local Group. 

Given the detection of dark microlensing objects, 
as well as the need for significant amounts of dark 
baryonic matter somewhere (Copi, Schramm, & Turner 



1995; Fields, Kainulainen, Olive, & Thomas 1996) 



halo white dwarfs are a very conservative candi- 
date (e^, Larson [l98^ ; Ryu, Olive, & Silk |1990 ; 
Silk |1993| ). This is particularly so since red and 



brown dwarfs are apparently excluded as halo can- 
didates (BahcaU, Flynn, Gould, & Kirhakos |1994 



Graff fc Freese |1996| ). Here we make a specific though 
schematic model, and assess the plausibility of the 
white dwarf hypothesis.^ As we will see, the model 
can be made to work but not without some assump- 
tions (e.g., one must alter the halo initial mass func- 
tion). In any case, the model is eminently testable, 
and perhaps has already been tested by X-ray ob- 
servations in other groups. Indeed, should one find 
this model and the halo white dwarf hypothesis un- 
tenable, then it follows that the dark baryons and the 
MACHOs must take an even more exotic form. 

2. Local Group Properties 

We take the Local Group to have a mass ~ (3 — 
5) X W'^Mq (Fich & Tremaine |l99l| ). Its luminous 
component is dominated by two galaxies, one of which 



^Others have suggested that the microlensing objects might 
be remnants from an early Pop ulation III; see Fujimoto, 
Sugiyama, Iben, & HoUowell (1995). 
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(the Galaxy) has a visible mass ~ 7 x lO^^M©. The 
total mass of the Galaxy is uncertain and depends on 
the radius, i.e., Mtot ^ 5 x lO^-^M0(i?haio/5O kpc). 
This implies that the dark halo has a mass within 
50 kpc of '--^ 4 X IO^^Mq, some or all of which will 
be in the form of microlensing objects and other dark 
baryons. We consider models with up to 90% the total 
mass in nonbaryonic dark matter, at the group scale. 
The baryonic fraction at the halo scale, however, can 
be much less. 

In many galaxy groups, hot intergalactic gas is 
found (e.g., Mul chaey et al. |1993| , |1996| ; Pildis, Breg- 
man, & Evrard 1995; Ponman, Bourner, Ebeling, & 



Bohringer 199t ). This gas is a significant and some- 
times dominant component of the baryonic mass. The 
hot gas-to-galaxy mass ratio ranges from 0.3 — 3. The 
gas also contains metals, with Z ~ (0.1 — 0.2)^0. 
While the ROSAT metallicity determinations are un- 
certain (Davis, Mulchaey, Mushotzky, & Burstein 



1996 ), recent measurements with ASCA (Fukazawa et 
al. 1996| ) suggest that the gas metallicities in groups 



could have a larger spread than those of clusters. In 
any case, the metallicity is clearly not primordial, in- 
dicting that a significant fraction of material has been 
processed in stars before they are incorporated into 
the intragroup medium. 

If the gas is in hydrostatic equilibrium with the 
gravitational potential of the group, then the tem- 
perature distribution determines the total mass. Fur- 
thermore, the observed temperatures tend to tightly 
cluster around T = 1 keV for most groups observed 



thus far (Mulchaey et al. 1996). This implies that the 
total masses are similar, with most observed group 
masses lying in the range (1.5 — 2.5) x 10^'^ Mq. This 
is a tighter range than the span of luminous mass in 
galaxies and gas. 

However, it is not clear that all optically identified 
groups evidence hot gas. Several analyses of ROSAT 
observations have found that the presence of de- 
tectable gas is strongly correlated with the morpholo- 
gies of the constituent galaxies (e.g., Pildis, Bregman, 
& Evrard 1995). The trend is very similar to that 



found in galaxy clusters: the hot gas mass is cor- 
related with the presence of early- type (E and SO) 
galaxies. Indeed, Mulchaey et al. ( 1996 ) emphasize 
that there is at least one bright (L_b > 5 x IO^^Lq) 
elliptical galaxy in every group for which hot gas has 
been detected. 

On the other hand, Ponman, Bourner, Ebeling, & 
Bohringer ( |1996| ), also using ROSAT data, have re- 



cently claimed the positive detection of hot gas in 
spiral-dominated Hickson compact groups. They find 
that the gas in these groups has a lower tempera- 
ture (~ 0.3 keV), and thus, a lower surface brightness 
than that in groups with early-type galaxies. The dis- 
crepancies between these results and those of previous 
groups attests to the difficulty in trying to measure 
or put limits on such a relatively dim diffuse com- 
ponent. Ponman et al. (1996) note that differences 



between their results and others trace to details of 
the analysis, e.g., subtraction of galactic and back- 
ground emis sion. As Davis, Mulchaey, Mushotzky, 
& Burstein (1996) point out, these are not always 
straightforward issues. Hence, we regard the issue of 
hot gas in spiral-dominated galaxies as presently am- 
biguous. 

In the Local Group specifically, there has not been 
direct obse rvatio n of hot intergalactic gas. However, 
Suto et al. (1996) have argued that such a component 
is allowed within current direct limits. Indeed, they 
suggest that this is a source of the excess low-energy 
component in the diffuse X-ray background. Suto et 
al. ( |1996D model a gas distribution that could lead to 
the excess radiation; their distribution implies a total 
mass in hot intergalactic gas of about ~ 3.5 x IO^^Mq. 
Of course, the existence of hot gas in spiral-dominated 
groups is a central premise in this scenario. Pildis & 
McGaugh (1996) show that the observational limits 
on such gas (if it is similar to other spiral-rich groups) 
require any Local Group gas to have too low a mass 
to provide the soft X-ray background. 

3. Hierarchical Collapse Model 

Motivated by hierarchical clustering scenarios for 
structure formation, we compute the evolution of a 
hierarchy of self-similar mass scales. The three scales 
are (1) protogalactic clouds, all of which reside in (2) 
galaxy halos, themselves moving within the (3) group. 
Within a spiral protogalaxy, the clouds merge to be- 
come ultimately the disk and bulge. However, we do 
not distinguish the disk and bulge. Their formation 
does not significantly affect the halo or group evolu- 
tion once the last merging has occurred and the star 
formation (and gas outflow) in the halo has dimin- 
ished. 

The dynamics of the three components are de- 



scribed (Mathews & Schramm 1993) as the radial evo- 
lution R{t) of a spherical overdensity, from its initial 
expansion (starting at t = 10* yr) through the de- 
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parture from Hubble flow and collapse to a fixed final 
radius. The halo collapse is halted at 50 kpc. 

For each component in the structure hierarchy, we 
follow the evolution of matter in the form of gas, stars 
and remnants, and possibly also non-baryonic dark 
matter. We also compute the helium and metal evo- 
lution, and follow the temperature of the hot gas. Our 
model is schematic but contains the necessary features 
for testing and constraining our basic hypotheses. Its 
structure (summarized pictorially in figure ^) is as 
follows. 

3.1. Level 1: Protogalactic Clouds 

In this simple model we assemble galaxies from 
a distribution of protogalactic fragments within ex- 
panding and contracting galactic halos. At the level 
of the protogalactic clouds, we consider the system to 
be composed of three components: cold star-forming 
molecular clouds; heated ejecta from the clouds; and 
stars and remnants. We assume that non-baryonic 
dark matter (if there is any) is unimportant at this 
scale. 

Let us first consider the cold star-forming gas and 
hot gas components. These components are most 
dramatically affected by star formation and merg- 
ing. Stars form from the cold gas and return a frac- 
tion of their material as hot ejecta. This hot stellar 
ejecta further mixes with and heats the local cold gas 
mass into the hot component. In addition, mate rial 
is heated during mergers (e.g.. White et al. 1993 ) as 
the relative kinetic energy of the merging clouds is 
converted into internal energy of the merged system. 
Also, cooling of the hot gas component returns mate- 
rial back to cold star forming regions. 

We write the coupled equations for the evolution 
of the cold and hot gas components for an average 
cloud experiencing all of these processes as: 
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The first term on the right hand side of equations (||) 
and (0) describes the formation of stars and ejecta. 



TZ is the usual (e.g., Tinsley 1980| ) returned fraction 
of material from stars in the instantaneous recycling 
approximation. It depends on the initial mass func- 
tion. The hot ejecta is assumed to sweep up and mix 
with the cold cloud material; the resulting mixture 
contributes to the hot gas component. We therefore 
include an additional factor of ?7mix describes the num- 
ber of equal masses of local interstellar material that 
is heated into the hot component along with the stel- 
lar ejecta. The rate of incorporation of cold gas into 
new stars is ip{t). 

The second term on the right hand side in equa- 
tions (|^) and (^ describes the effects of mergers. 
Here, Amorgc is the merger rate per cloud. On average, 
during each merger, the cloud mass will double. The 
third term describes the cooling from hot to cold gas, 
with the cooling rate per particle given by Acooi . The 
last terms in Eq. (||) accounts for infall and outfiow of 
gas. The infall term min gives the rate at which hot 
halo gas is incorporated int to the hot gas component 
of the clouds. The outflow rate TOquJi we attribute to 
the evaporation of hot gas. The physics behind each 
of these terms is discussed in §3.4. 



For our purposes, we take the cold star forming gas 
to be truly cold gas with T < 100 K. The hot gas com- 
ponent is approximated by assuming that all gas at 
higher temperatures is isothermal and homogeneous 
within each of the levels of structure. 

If we combine equations (|l|) and (|^) we can recover 
the familiar instantaneous recycling equation for one 
gas component; 



^cold 



hot 



-(l-7^)^(^) + <f -mout , (3) 



where the effective infall rate, due to mergers as 
well as the influx of halo gas, is TO?f = {m'^'^ + 
"T'!?i*)Amcrgc + "^ini- The outflow rate is just mout = 

•hot 
'"■outl- 

As in MS93 we assume that protogalactic mergers 
disperse the stars and remnants into the halo. The 
evolution equation for the mass rriri in stars and rem- 
nants remaining in a cloud becomes: 



TOj-i = (1 - Tl)lpit) + (1 - K)Amorgc"^rl 



(4) 



Here k is a measure of the efficiency for the mergers to 
disperse stars and remnants. Specifically, the average 
fraction / of stars and remnants born at time to which 
survive merging up to time t is 



f = e 



I r* dt' 



merge 



(5) 
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We will take k = 1. 

Since stars return all of their ejecta into the hot 
component, the metallicity in the cold star forming 
gas is only indirectly enriched by cooling from the 
hot component. Thus, from Eq. we write for the 
evolution of the total mass in metals in the cold star 
forming gas (Zf 



,hot 
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This reduces to a simple evolution equation for the 
metallicity of the cold gas, 



^cold ^^hot 



7Cold\ 



A, 



cool " 



^hot 



^cold 



(7) 



Clearly, the equilibrium metallicity of the cold star- 
forming component is equal to the hot component 
metallicity. This equilibrium is, however, only achieved 
after a time given by the cooling time scale, Acooi^^ 
times the ratio of cold to hot gas masses. 

By an analogous process to the derivation of Eq. 
(j^), the evolution of metallicity in the hot gas com- 
ponent can be written 



i'hot 



[VZ + (1 + ?/mix)7^(Z^°'^ - z\°% 
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where yz is the mass fraction of newly synthesized 
material in the ejecta. As usual, yz depends on both 
the initial mass function, and on the stellar nucleosyn- 
thesis yields as a function of mass (and metallicity). 

By a similar derivation, the temperature of hot 
gas in the halo is determined from an energy bal- 
ance equation. This leads to an evolution of internal 
energy per unit mass ei in the hot gas. 
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^We explicitly show only the metal evolution here and below, 
but from this the expressions for the helium evolution follows 
trivially. 
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where egN in the average energy per unit mass in 
all stellar ejecta averaged over an appropriate initial 
mass function. This term is dominated by supernovae 
whose energy release is E'oj ^ 10^^ erg per supernova. 
Averaging over the initial mass function 0(to) then 
gives 

J dm (f>{m) Ef.j{m) 



(10) 



J dm (j){m) TTioj 

which is independent of the normalization convention 
of (p. For a typical initial mass function, this gives 
values of order esN ~ lO'^^ erg M~\ 

The quantity e°"' is the average energy per unit 
mass for the material exiting the cloud. For an ideal 
gas, the temperature is simply related to ei by, Ti = 
(2ei/i/3fc7Vyi), where is the mean molecular weight 
of the gas and Na is Avagadro's number. 

The last term in Eq. (^) accounts for the p dV 
work done as the cloud expands or contracts; the ra- 
dial dependence is given by the collapse of a spherical 
overdensity. The cloud radius i?i is the average tidal 
radius 

1/3 

i?2 , (11) 
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where toi is the total average cloud mass 
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and 7712 is the average halo mass defined below (Eq. 
p^ ). We assume that the relative momentum of the 
merger goes into heating the cloud gas. Thus, we have 

— -Emerge /-^^^merge 



Emerge 



^?el/2. 



3.2. Level 2: Galactic Halos 

For our purposes, galactic halos are treated as a ho- 
mogeneous assembly of protogalactic clouds and hot 
gas, possibly having a component of nonbaryonic dark 
matter. The total mass of the galactic halo is then 



^2 = 7lcTOl + TO<r2 + TO,.2 + m 



NB 



(13) 



where 77c is the number of protogalactic clouds, 777g2 
is the mass of hot gas which has exited the clouds to 



5 



reside in the halo, and 771^2 is the mass of stars and 
remnants which have been dispersed from clouds into 
the halo. The mass of the nonbaryonic component, if 
present, is m^^. 

The equation describing the evolution of the (hot) 
gas component in the halo is: 



"7g2 = f^c ("^outl - Willi) 

— "^out2 + "^in2 J 



(14) 



where mout2 is the rate at which the hot gas is ejected 
from the halos as described below and rhins is the 
possible inflow of gas from the intragroup medium. 

The rate at which stars and remnants are injected 
into the halo from mergers can be inferred from Eq. 



merge 



(15) 



The equation governing the evolution of metallicity 
in the halos will be 



Z2 = (^r-^2; 

+ {Z3 - Z2) 
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(16) 



m„2 



where ric is the number of cold protogalactic clouds in 
the halo. Similarly the internal energy per unit mass, 
€2, of gas in the halo is determined from an energy 
balance equation: 
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The number of clouds decreases exponentially with 
the number of mergers: 



= -A 



merge '^c 



(18) 



3.3. Level 3: Hot Intra-Group Medium 



The equations governing the evolution of the intra- 
group medium are analogous to those for the galactic 
halos, but at the level of the group nonbaryonic dark 



matter may be a dominant contributor. Thus, for the 
total group mass we write, 



ma = nhTO2 + "T-gS + "t-. 
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(19) 



where rih is the number of galactic halos in the group, 
rrigs is the mass of the hot X-ray intragroup medium, 
and m^^ is the contribution from nonbaryonic dark 
matter. 

The evolution equation for the intragroup medium 
is then, 

rn^3 = n\,[m2out - TOin2) - mzout , (20) 

and msout is the rate at which the gas is lost from the 
group. 

The equation governing the evolution of metallicity 
in the intragroup medium is then 



Z3 — {Z2 — Z3 



mg3 



(21) 



and the energy balance equation is 

. /out ^ ^ nhmout2 

£3 = (£2 -£3; 

mg3 

- 1^3 - £3) 



(22) 



mg3 

— '\cool£3 — ^7^£3 I 

from which temperature of the intragroup medium 
can be inferred. 

3.4. Cooling, Star Formation, Mergers, and 
Mass Loss 



The evolution is given by equations ([i|-p2|). What 
remains is to specify the various input quantities, i.e., 
the rates for cooling, star formation, merging, and 
mass loss. 

The metallicity-dependent cooling rate Acooi is de- 
rived from the calculations of Bohringer & Hensler 
( 1989 ). For the temperatures appropriate for the 
hot gas, the cooling is dominated by brehmsstrahlung 
emission from electrons. These losses are written in 
terms of the cooling function A, which gives the en- 
ergy loss rate per unit density. Given the number 
density at a given scale Ui, one may then compute 
the energy loss rate per particle: Ei — n^A, and then 
the cooling rate Acooi = E/E. 

We allow star formation to be induced both by 
mergers and intrinsic star formation processes within 
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the clouds (MS93). In MS93 it was shown that m- 
trinsic star formation is more important during the 
subsequent evolution as material settles into a disk. 
To describe the merger-induced star formation for all 
of the merging substructures within a collapsing halo, 
MS93 proposed a schematic model of colliding viri- 
alized protogalactic clouds within an expanding and 
contracting halo. We adopt that formulation here. 

Specifically, we presume that the fragments virial- 
ize as they form. The collision rate per cloud within 
this virialized velocity distribution can be written. 



An 



K - \)avlV 



(23) 



where n is the number of protogalactic clouds within a 
volume y , a is an average collision cross section, and v 
is the virial velocity {y ~ [0.4GMh/i?h]^/^), where Mh 
is the total gravitational mass of a galactic halo, and 
it!h is the radius of the galactic halos, approximated as 
a collapsing spherical overdensities. We define halos 
as those regions which evolve to become independent 
gravitationally bound ensembles of gas and stars at 
the present time. All halos themselves are viewed 
as the result of merging internal structure. Hence 
we only describe merging within the halos and not 
merging between halos. This is largely a matter of 
semantics. 

The number of protogalactic clouds decreases ex- 
ponentially with the integral of the merger rate (Eq. 



= nc(0) e'/o •^'°°>^g<="^ , (24) 

where the initial number of clouds, 7ic(0), is given by 
the ratio of the total initial halo baryonic mass to 
the initial cloud mass. Thus, nc(0) ~ 10® for this 
schematic model. For the merger cross section a we 
use cr = 7ri?^, with a tidal radius i?t = 2i?i to include 
gravitational effects (Binney & Tremaine 1987). 



Near the end of the collapse of the halos, the com- 
bined effects of conservation of angular momentum 
and heating will dissipate the radial motion. In our 
schematic model, we approximate these effects by 
halting the halo collapse at a size of a typical present 
dark-matter halo, (i?2 ^ 50kpc). The cloud radii are 
fixed by Eq. (|ll|), and the group collapse is halted 
when the radius reaches the present size of the Local 
Group, i?3 ~ 700 kpc. 

The merger-induced stellar birth rate thus is taken 
as proportional to the mass of gas mass participating 
in mergers per unit time. The intrinsic, quiescent 



star formation rate is taken as oc p" 



)id' 



^cold 



71 = 1/2 in our case. The total star formation rate is 
the sum of the two terms: 



cold 



(25) 



Note that we have parameterized the strength of 
merger-induced star formations in terms of the di- 
mensionless efficiency amrg- Typical numbers for 
1% (MS93). The coefficient of the qui- 
escent piece is taken as /3q = /3qAo, where Aq = 
1.7 X 10-V(Mq kpc-^)-i/2 Gyr"^ Typical values 
of the dimensionless scaling /3q are 1 — 5 (MS93). 

Finally, we must specify the rate at which hot gas 
is ejected from, or falls into, the clouds, halos, and 
the group. We first calculate the outflow due to evap- 
orative mass loss. We assume that the hot gas is dis- 
tributed homogeneously within an object. The rate 
of loss of hot gas is then just given by the fraction of a 
Maxwellian thermal distribution of velocities in excess 
of the escape velocity at the surface of the structures 
at all three levels. Thus, we write. 



hot 



^ (^(> ^e)) hot 
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(26) 



Where, i denotes clouds, halos, or the group, and 
Cioss < 1 is a dimensionless scale factor depending 
on the geometry of the cloud; Cioss = 1 for a sphere. 
The factor (t;(> Wosc)) is the average velocity of all 
particles above the escape velocity v\^^ = 2GMi/Ri: 



{v{> Vcsc)) = 



/Cc hiB{v) 
dv fMB{v) 

{l + x,)e~^'VT (27) 



where /m b is the Maxwellian distribution. The aver- 
age thermal velocity is = 2kTi/ fimp, and /xm^ is 
the average mass of a gas particle. Xi — GMi/ RikTi 
is the ratio of gravitational binding energy to thermal 
energy, and is small typically. 

Similarly, the average energy loss per unit ejected 
mass of material is 



{ev{> Vcsc)) 

{v{> Vcsc)) 

1 + Xi+ xf /2 2kTi 

1 + Xi /ITOp 



(28) 



where 



which for small Xi gives e°"* ~ 2kTi/finip = 46^/3. 
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Using similar derivations, we have computed in- 
flow of halo (group) gas into the clouds (halos). We 
note that this inflow has two components. First, some 
material is accreted simply due to the motion of the 
halo or cloud through the surrounding medium-the 
finite size of these substructures will lead to the cap- 
ture gas that falls within their geometric cross section 
a = TrRf. Second, the hot gas particles surround- 
ing the substructures will accrete due to their ran- 
dom thermal motion. We assume that these particles 
quickly equilibrate with the ambient gas. Together 
these two effects lead to an infall term of the form 



Clos: 



Ri 



i+i 



hot ( T 



R 



(29) 



with v1 = Gm. 
locity. 



/ Ri+i the cloud (halo) circular ve- 



3.5. Initial Mass Function and Stellar Yields 

The initial mass function (IMF) is important for 
our scenario in several respects. As we have noted, 
it affects the parameters of our model, namely the 
returned fraction TZ, yields j/y and yz, and specific 
energy injection esN- Moreover, the crucial impact of 
the IMF for our scenario is that it determines the ratio 
of low mass stars to high mass stars in the halo. Low 
mass stars (m < O.OMq for halo metallicities) will 
still be burning today. While these stars are faint, 
they would be detectable, and heavily favored by an 
IMF similar to that of present disk stars. As a re- 
sult, such stars would today outnumber by far the 
halo remnants, and their net contribution to the halo 
luminosity would be large (Ryu, Olive, & Silk 1990 ). 
Thus, remnants cannot be a significant component of 
the halo if the halo IMF is similar to that inferred 
for disk stars. To allow for a significant population 
of halo white dwarfs, they must have been formed 
from an IMF which strongly favored the formation of 
intermediate- to high-mass stars over low-mass stars 
(Ryu, Olive, & Silk |199C| ). 

Given the different physical conditions during the 
formation of the halo (e.g., frequent mergers, higher 
temperatures, lower metallicities), it is at least plau- 
sible that the star forming process then was different 
than that of disk stars. Indeed, there are arguments 
for an early IMF that is skewed towards higher masses 
(see Silk |l993| ; Bond, Arnett, & Carr |1984| ; Adams & 
Fatuzzo |1996|). 



vational (Miller & Scalo 1979|) a nd theoretical ar- 
guments (Adams & Fatuzzo |1996| and refs therein): 
ln0(lnTO) = Inc^o ~ [ln^(m/TOc)/2(T^]. We investigate 
the effect of variations of the centroid mass rric and the 
dimensionless width a. We take the IMF bounds to be 
(0.1,100) Mq, w ith a black hole cutoff at 18MQ(Brown 
& Bethe p^ ). In fact, our results for metal yields 
are insensitive to these limits. With these parameters, 
we obtain a returned fraction TZ = 0.375. 



Therefore, we have chosen the simple log-normal 
IMF parameterization, as suggested by both obser- 



The stellar yields are taken from Maeder (1992). 
We follow both helium and metallicity Z. These 
abundances are particularly powerful constraints when 
used together. The helium yields are mainly from in- 
termediate mass stars (to < 8Mq), while the metal 
yields come primarily from high mass stars. Thus, 
the IMF must strike a balance between the two mass 
ranges to avoid an inappropriate ratio of helium to 
metals. With our IMF, we find yy = 0.0139 and 
yz = 0.0354. 

4. Results 

In this hierarchical merging picture we take the 
clouds to be the initial building blocks of all structure. 
Thus, while we assume the clouds to be comprised 
initially of hot and cold gas, we assume baryons in 
the halos and the Group to be initially within clouds 
only. Thus, by definition the initial gas in the halos 
and group is zero. The initial metallicity is zero, while 
we take a primordial helium abundance of y = 0.235. 

We have run the model for plausible ranges of the 
input parameters. The two models we present here 
are: (1) a "best-fit" model which has appropriate fi- 
nal masses and metallicities , while a lso optimally sat- 
isfying other constraints (§§4.4-4.6); and (2) a "max- 
imum remnant" model which has the highest possible 
halo remnant mass without violating the constraints. 
Parameters of these models are summarized in Table 
1. Unless otherwise noted, numerical results will be 
given for the best-fit model. 

For all models, the qualitative results are as fol- 
lows. There is generally a high initial merger rate, 
which lasts for < 1 Gyr. This reduces the number 
of protogalactic clouds from 10^ to about 100 while 
producing a burst of star (and hot gas) formation. 
The star formation rate and cloud number evolution 
appear in Figure ^j. Some of the stars are ejected into 
the halo by the mergers, and on a longer timescale, 
a wind is ejected — first from clouds, then from halo, 
and ultimately from the Group. When the halo col- 
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lapses, after 5 Gyr (Mathews & Schramm 1993 ), there 
is an additional burst of merging and star formation. 
The remaining clouds coalesce into what will eventu- 
ally become the galactic bulge and spiral arms. With 
the collapse there is also heating and a reinvigorated 
wind. Subsequently, all remaining (hot, metal-rich) 
halo gas is ejected into the intergalactic medium. The 
evolution of halo and group masses is given in Figure 



4.1. Disk and Bulge Formation 

The clouds begin with the Jeans mass at recom- 
bination, IO^Mq. They merge to form the proto- 
disk-l-bulge, with a mass 8.1 x 10^°Mq. It is encourag- 
ing that, at the end of the halo collapse, most material 
in the proto-disk is in cold gas which will form disk 
stars. The evolution of metallicities is given in Figure 
^ As seen in the figure, the metallicities grow rapidly 
in the initial burst, then remain fairly constant un- 
til the halo collapse at 5 Gyr. Then the metals rise 
again as the halo collapses and the disk and bulge are 
formed. Hence, the halo evolution provides an initial 
metallicity for the gas of the protodisk. This corre- 
sponds to an "initial enrichment" of the disk material, 
and so avoids overproduction of metal poor disk stars 
(i.e., the disk G-dwarf problem). 

4.2. Halo Formation 

In the best-fit model, portrayed in Fig. the 
galaxy (=clouds-|-halo) begins with a mass of 1.35 x 
lO^^Mo, of which 1.1 x IO^^Mq is baryonic (all 
in the clouds). Even higher fractions are allowed, 
as discussed below (§||). The final galaxy mass is 
5.0 X lO"Af0, of which 2.5 x lO^M© is baryonic; thus 
50% of the total Galactic mass is in baryons. Some of 
these baryons constitute the disk-|-bulge, as described 
above. The baryonic mass of the halo remnants alone 
(excluding the disk-|-bulge) is 1.7 x lO^^M©, and the 
halo nonbaryonic mass accounts for 2.5 x 10 ^"'^M©. 
Hence, 40% of the dark halo is in remnants, consis- 
tent with the microlensing observations. 

When combined with mass loss from the other 
galaxy in the Local Group, there is a total loss of 
1.7 x IO^^Mq of gas from both halos into the intra- 
group medium. The galactic wind is thus very ef- 
ficient, as it must be to remove the ejecta that ac- 
companies the remnant production. The remnants 
themselves are mostly white dwarfs, with 12% neu- 
tron stars. 



4.3. Group Evolution 

In the best-fit model, the Local Group begins (Fig. 
PD)with a mass of 5.6 x lO^^M©, and a baryonic mass 
of 2.2 X lO^^M©. At the end of the simulation the 
group has a total mass 4.3 x lO^^M©, 17% of which is 
baryonic. Of the baryonic group mass, 2.9 x lO^^M© 
resides in hot gas. The intragroup gas temperature 
is 0.26 keV, just at the limit of ROSAT sensitivity. 
The group as a whole loses 1.4 x lO^^M© of gas to 
the intergalactic medium; thus about 64% of the ini- 
tial baryonic mass in the group is ejected later into 
intergalactic space. 

It is remarkable that the remnants in the halos, 
turn out to represent a small fraction of the initial 
baryonic matter. The galactic winds required to re- 
move the stellar ejecta co-produced with the white 
dwarfs prove strong enough to remove most gas from 
the group itself. This amount of hot (ionized) ma- 
terial is consistent with Gunn-Peterson limits on the 



intergalactic medium if it does not cool further (§5.1). 

Thus, in the best-fit model we find the dark halo 
to be 40% microlensing objects. These take the form 
of stellar remnants, 88% of which are white dwarfs, 
and the rest neutron stars (and perhaps black holes). 
In this model the copious production of hot intra- 
group and intergalactic gas is a natural consequence of 
white dwarf-dominated halos. We produce a present 
mass of intragroup gas of 2.9 x lO^^M©. This cor- 
responds to about 37% of the baryonic mass of the 
Local Group. This mass is encouragingly near the 



value (3.5 x 10"Mq) imphed by the Suto et al. ( |l996D 
model for the excess diffuse X-ray background (§y). 
(But recall the observational controversy §|[) 

The X-ray luminosity of the halo is 1.5x10'*'^ erg s~^. 
Let us assume for simplicity that the earth is at the 
center of this emission, and that the gas extends ho- 
mogeneously to the edge of the Local group (i?3 = 700 
kpc in this model). This leads to a total diffuse back- 
ground of about 4 X 10~^ erg s~^ cm^^ sr^^. This is, 
at its peak energy, about 3 orders of magnitude be- 
low the observed diffuse soft X-ray background. We 
therefore find the diffuse Local Group radiation to be 
negligible. The contribution is thus much lower than 
that suggested by Suto et al. (1996), who postulate a 
similar (but larger) intragroup gas mass, but posit a 
temperature of 1 keV as opposed to 0.26 keV in this 
model. 
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4.4. Mass budget 

The most basic constraint on our model is the re- 
quirement that it reproduce the observed mass and 
metal budget of the Local Group. The final masses 
determine the initial masses, but the metallicity val- 
ues constrain how to get these. Consider the to- 
tal metal mass production in our model. The to- 
tal initial baryonic mass is 2.2 x 10^^. Of this, the 
mass processed into stars and remnants is Mrem = 
4.95 X IO^^Mq. The total mass into ejecta is just 
7^/(l - n)/Mrem = 3.0 X IO^Mq. The total metal 
mass produced is yzMe-^ = 1.1 x lO^^M©, and the 
total new helium mass is uyM^i — 2.8 x IQ^^Mq. 

Now, the average new metals from a given star rep- 
resents a high fraction of that star's ejecta: yz = 
0.014 ~ 2/3 Zq. However, this metal- rich ejecta is 
diluted in several ways. First, it is mixed with 77 
times its mass as it leaves the cold star-forming re- 
gions. Then it is mixed into the hot cloud gas, which 
quickly evaporates into the halo and on to the intra- 
group medium. Much gas is even evaporated from the 
Local Group itself, becoming part of an intergalactic 
(extragroup) medium. Thus, the final mean metallic- 
ity must be averaged over all of the baryons (including 
the ejecta from the Local Groups). One finds that the 
global average mass fraction is just 0.0049 ~ 0.25^©, 
an acceptably small value. 

The key point here is that star formation does not 
occur in all regions containing gas (as is often as- 
sumed in simple chemical evolution models) but only 
cold molecular cloud cores. Consequently, the stars 
represent a small fraction of the total baryonic mass. 
Furthermore, whereas the stars remain in the galax- 
ies, the gas is dispersed at all levels. The halo can 
thus contain many remnants but not a large amount 
of metals. 

4.5. Nucleosynthesis 

Nucleosynthesis provides an important additional 
constraint on t his sc enario (Ryu, O live, fc Silk 199C| ; 
Chariot & Silk |1995| ; Hegyi & Olive |1986| ). The burst 
of star formation which produces the white dwarfs 
must not overproduce metals and helium. As noted 
above, the metallicities at the different scales are rea- 
sonable. Furthermore, we may test not only the aver- 
age metallicities, but also the halo metallicity distri- 
bution. Figure || compares our calculated halo stel- 
lar metallicity distribution with the observed globu- 
lar cluster distribution. The agreement is good for 



the most metal-poor {[Z] < —1), and thus the oldest, 
portion of the population. 

The halo metallicity distribution strongly constrains 
the total halo remnant mass in our model. The posi- 
tion of the peak constrains the net star formation rate 
amplitude. If the star formation rate is too large, the 
peak is shifted too high. Thus, the need to reproduce 
the peak at [Z] ~ —1.5 limits the mass processed 
into stars, and so leaves fewer stars available for the 
halo. From the of a fit to the halo metallicity 
distribution we estimate a 2a upper limit of a 77% 
remnant contribution to the halo. Consequently, it 
seems likely that the halo contains at least some non- 
baryonic dark matter. In addition, the shape of the 
distribution is largely sensitive to the relative star for- 
mation rate, and to the degree of mixing between the 
stellar ejecta with the cold gas. 

The IMF shape is of course important in deter- 
mining the metallicities. However, contrary to Ryu, 
Olive, & Silk, we find no significant constraint on the 
IMF upper limit, essentially because the ejection of 
gas removes some of the metals produced by high- 
mass stars. This is even more so if there is a cutoff 
for black hole formation as low as 18 Mq (Bethe & 



Brown 1995). This is 



4.6. Halo Luminosity 

Several sources of luminosity also provide impor- 
tant constraints. 

(1) Low mass stars in the halo are long lived and 
can lead to an unacceptab ly hig h halo mass-to-light 
ratio (Ryu, Olive & S ilk 199C| ; Richstone, Gould, 



Guhathakurta, & Flynn 1992); however, in our model, 
the mass to light ratio is comfortably above these au- 
thors limit M/Lb > 500 for i?haio = 50 kpc. We 
determine the luminosity contributed by all stars still 
burning in the halo by integration of stellar luminosi- 
ties, the past star formation rate, and initial mass 
function, without the instantaneous recycling approx- 
imation. We find a mass to total luminosity ratio of 
M/Ltot = 760, M/Lv = 1300, and M/Lb = 3500. 

(2) Light from the bursts of star formation can 
lead to a large diffuse background (Chariot & Silk 
1995| ; Zepf & Silk 1996). Using a population synthe- 
sis model we have preliminarily computed this back- 
ground assuming the starlight to be unscattered after 
its emission. These will be discussed in Mathews et al. 
(1996), but results are consistent with observed con- 



straints. We exceed the Chariot & Silk (1995) limit 
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of a 10% remnant halo because most of the elements 
are formed at very large redshift (z ~ 10). 

(3) The halo white dwarfs will contribute to the 
white dwarf luminosity function (e.g., Adams & Laugh- 
lin |1996| ; Chabrier, Segretain, & Mera |1996D . There 
is a new constraint from Flynn, Gould, & Bahcall 
(1996), which requires halo white dwarfs to have an 
luminosity My > 18.4. for y - / > 1.8. 

The white dwarf luminosity function in our model 
is given in Figure ||. The two peaks corresponding to 
the two bursts of star formation. The strongest con- 
straint comes from the more recent peak produced at 
the epoch of halo collapse. These white dwarfs have 
cooled to present luminosities of ~ 1O~^-^L0. This 
would correspond to a bolometric magnitude of 18.4. 
However, the bolometric correction (Liebert, Dahn, 
& Monet 1988| ) , would probably reduce the visible lu- 
minosity to Aly > 20, well below the Flynn, Gould, 
fc Ba hcall limit (similar to the findings of Kawaler 
1996 ). The low luminosity of the white dwarfs de- 
rives in part from their very old age, but also because 
the population comes from a shifted IMF. Since the 
progenitors are more massive stars, they died sooner 
and so the remnants have had longer to cool. Indeed, 
our adopted IMF centroid of 2.3Mq is that of Adams 
& Laughlin ( 1996| ), who chose it specifically to obey 
the luminosity function constraints in the disk. 

4.7. Parameter Sensitivity 

Our results depend sensitively on several of the 
model parameters. Parameter values for our best-fit 
model are summarized in column 3 of Table 1. Among 
the most important parameters are those that deter- 
mine the shap e of the IMF. We find a good fit for a 
centroid (c.f. §3.5) rric = 2.3Mq, and a width (7 = 1. 
Note that this centroid is that of Adams & McLaugh- 



lin (1996). However, they advocated a tighter width 
to avoid metallicity problems which we do not require 
in our model because metals are efficiently ejected in 
the hot wind. Indeed, the width we use is essentially 
the present value (Miller & Scalo 1979| ). 

Other IMF parameters are allowed, but to obey the 
observational constraints, these must be near the ones 
we have adopted. As discussed above (§^.5]), the cen- 
troid must be > IMq to avoid significant production 
of long lived, low mass stars still burning in the halo. 
However, the centroid must not be too high to avoid 
untenably large metal yields. The width must also be 
not so large that it allows too many low mass stars. 



but not so narrow that metal yields are too small and 
helium yields too large. Thus, the need for a dark 
halo, and for a reasonable nucleosynthesis, drive the 
IMF parameters to the range of those we have chosen. 

The other important parameters involve star for- 
mation and ejecta. The star formation parameters 
(eq. ^ and Table 1) are reasonable and consistent 
with previous results (MS93), derived from requiring 
consistency among various comoschronometers. The 
quiescent term dominates, controlling the mass bud- 
get, but the merger term is also important. The ratio 
of the merger to quiescent contribution controls the 
position of the peak in the halo metallicity distribu- 
tion (Fig. |). 

For the ejecta, the most important parameter is the 
degree of mixing Tymix- A higher value leads to more 
dilution of the ejecta, which leads to lower metallicity 
and temperature in the hot gas. The observed lower 
temperatures and metallicity of the hot gas in groups 
demands that rj^^ix > 1. 

5. The Remnant Contribution to the Dark 
Halo Mass 

Given that (1) there is good evidence for dark mat- 
ter in the Galactic halo, and that (2) microlensing ob- 
jects are the first positively detected dark matter in 
the halo, a key question arises: how much of the dark 
halo is in MACHOs? In this context it is important to 
note that our best model does not give an all-remnant 
halo. Such a model is attractive for its simplicity, 
but we find that our best fit prefers only ^ 40% of 
the dark halo mass in remnants, with the balance 
in non-baryonic material. The strongest constraint 
on the white dwarf halo fraction is the halo metal- 
licity distribution. Increasing the white dwarf frac- 
tion requires increasing the star formation rate; this 
leads to h alo metallicities whose distribution peaks 
too high (§4.7). The "maximum remnant" model is 
that having the 2a upper limit for the halo fraction 
(as set by the halo globular cluster metallicity distri- 
bution). It's parameters are summarized in Table 1; 
in it the halo remnant mass fraction is 77%. This up- 
per limit can be raised somewhat if we increase the 
stellar merger dispersal efficiency k (see §3.2). While 
it is difficult to increase the halo remnant fraction, it 
is easy to decrease it. If the winds are less efficient 
(e.g., if Cioss < 1), then more material is recycled in 
the clouds. Thus a lower star formation amplitude is 
required and fewer stars are available to eject at early 
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Table 1: Model Parameters Adopted 



Model Feature Parameter Best-Fit Max Remnant 



IMF rric 2.3 Mq same 

£7 1.6 same 

Star Formation amrg 0.007 0.002 

/3q 1 2.5 

Ejecta r?inix 8 8 

Closs 1 1 



times. We find a minimum halo remnant fraction, for 
Closs < 0.08, of ^ 1%. 

Our findings are to be compared to recent analysis 
of the MACHO data (Alcock et al. |1996D . The mi- 
crolensing results are still very model-dependent, but 
it is provocative that they find a best fit for ^ 50% 
of the halo being made up of O.5to'2-^0 objects. The 
question of whether one can put stronger constraints 
on the MACHO halo fraction in our model leads to 
the issues raised in the next section. 

5.1. X-ray Groups: a Key Constraint 

As indicated in there is currently large uncer- 
tainty regarding X-ray observations of spiral domi- 
nated groups. However, while the observations are 
ambiguous, what seems more clear is that spiral-rich 
groups do not contain hot gas with T ^ 1 keV. That 
is, these systems do not evidence the same kind of 
diffuse emission as groups with early- type galaxies. If 
hot gas exists in these groups, it must have either low 
mass (< (1—3) x lO^^M©), or a low temperature, near 
or below the ROSAT threshold Tth < 0.3 keV. Indeed, 



Mulchaey, Mushotzky, Burstein, & Davis (1996) sug- 
gest that these systems do have cool (T < 0.3 keV) 
gas, and argue that this may have been detected as 
high-ionization quasar absorption lines. Furthermore, 
given the morphological similarity of the bulge and 
halo with early-type galaxies, it would be surprising 
if no hot gas were found in these systems. 

At any rate, the observations seem to rule out that 
the gas is both massive and hot in spiral- rich groups. 
However, in our model the temperature and the gas 
mass are related, as follows. 

After the processing of one stellar generation, the 
ejecta to remnant ratio is Tl/{1 — TV). Mixing with 
unprocessed cold material gives a total hot gas to 
remnant ratio (1 -I- J7mix)^/(1 — TV)- In our model, 
n = 0.38, so Mgas > Mrcm for 77mix > 1-7 (which 
holds in our case) . This already means that the mass 



of ejected gas must be larger than the halo baryon 
mass (assuming that most hot gas escapes). Con- 
sequently, if the limits on gas in spiral-dominated 
groups apply-i.e., if the gas were to have T > 0.3 
keV-then the mass in halo white dwarfs must be less 
than the X-ray limits on the gas mass, ~ lO^^M©. 
This would be only a small component of the halo. 

To allow for a significant gas mass, therefore, one 
must demand that it avoid the ROSAT limits because 
it is cool: T < 0.3 keV. However, the gas temperature 
is related to the gas-to-remnant ratio, as we now show. 
The total energy in the ejecta is 7?.esN-Mrem- This 
implies that the energy per unit mass of ejected gas 

is egas = 7^esNA^rom/^4as = eSN(l " 7^)/(l -|- 77,„ix)- 

Finally, using the scaling Tg^ = (2/3) /impfgN ~ 2 
keV, we have Tgas = 2 keV(l - 7^)/(l + ?7mix)- That 
is, the gas temperature diminishes with 77inix- But we 
already have seen that having a significant fraction of 
the halo in remnants, we need the concomitant halo 
gas to be cool. Using the Ponman et al. ( |l996|) value 
Tgas — 0.3 keV, and with TZ — 0.39, one finds ?7inix — 
3. In our detailed model which includes reprocessing, 
we find a larger value, 77 = 8. Thus, cooler gas results 
if '7inix > 1 , but only if the bulk of the baryons are in 
the gas. 

5.2. Cosmological Baryon Fraction 

In the preceding section we see that our best-fit 
model finds that a large fraction of baryons, initially 
in the group, are ejected as an intergalactic medium. 
This can be reconciled with the total cosmological 
baryonic budget, as follows. In our scenario, the cos- 
mic baryonic inventory is Sl^ = r^disk+buigc + f^hak> + 
Kg + ^!gm ^ KIZ + Kg + ^;^fGM (writing the uni- 
versal density in component i in units of the critical 
density: fi^ = p.j/pcrit, where pcrit = SHq/SttG). If 
most baryons reside in the intragroup and intergalac- 
tic diffuse gas, we therefore require ri^c + ^igm — 
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Observationally, galactic rotation curves imply that 
galaxy halos have fihaio ~ 0.02/i(i?haio/50 kpc), where 
-Rhaio is the (unknown) radius of the dark halo (e.g., 
Peebles 1993| ). On the other hand, primordial nucle- 
osynthesis calculations give Hb — (0.015 ± 0.005)/i~^ 
(Copi, Schramm, & Turner |1995| ; Fields et al. |1996D . 
Taking the ratio, we have Jlhaio/^^B ^ 1.3h^. Conse- 
quently, for h < 0.9, r^haio < ^B', that is, for reason- 
able values of h, some baryons are likely to be non- 
galactic. Indeed, for a low Hubble constant {h = 0.5), 
the bulk of the baryons are nongalactic. Furthermore, 
these relations assume that the halos are entirely 
baryonic {not the case in our best model). Any non- 
baryonic component in galaxy halos only strengthens 
the argument. 

Thus, we see that the X-ray observations of groups 
can provide strong constraints on remnants in the 
halo. Of course, this assumes that the Local Group is 
like other poor groups. If so, and if the X-ray data is 
reliable, then these data may provide a key constraint. 
As we have noted, if there is only a small X-ray gas 
mass in the Local Group, then there can be few rem- 
nants in the halo. Taking the more optimistic view, if 
the Ponman et al. ( |1996| ) result is correct, then there 
could be a large amount of (as yet unobserved) hot 
gas in the Local Group, as required in our model. 

6. Conclusions 

We have shown that, without violating constraints 
posed by luminosity and nucleosynthesis considera- 
tions, one may construct a plausible model in which 
the dark halo of the Galaxy contains a significant 
fraction of white dwarfs. These may have already 
been detected in halo microlensing events towards the 
LMC, and might also be detected via their luminos- 
ity function. The same bursts of star formation which 
produced the white dwarfs also led to hot, metal-rich 
intergalactic gas, some of which may still reside in the 
Local Group. This hot gas could be detectable via its 
X-rays, and by distortions in the cosmic microwave 
background radiation (Suto et al. 1996). 

Thus, the predictions of the model are testable. If 
the halo is comprised of white dwarfs then there must 
be a background of hot. X-ray emitting gas in the 
Local Group. Conversely, if there is metal-rich hot 
gas in the Local Group, then a significant fraction of 
the halo mass must be in remnants. Clearly, further 
searches for both of these are warranted. 

Furthermore, if our galaxy formation scheme is in- 



deed universal, then hot gas production and ejection 
should be a ubiquitous aspect of halo formation. Con- 
sequently, X-ray observations of other systems could 
provide a key constraint on our model. In partic- 
ular, our model can be directly tested by observa- 
tions which can unambiguously confirm or deny the 
presence of hot gas in other spiral-dominated groups. 
Also, if white dwarfs are ubiquitous in galactic ha- 
los, then they may lead to detectable infrared pro- 
files in edge-on galaxies, which may already have 
been observed (Barnaby & Thronson 1994; Sackett, 
Morrison, Harding, & Boroson 1994; Leq ueux, Fort, 

|1996| ; Lehnert & 



Dantel-Fort, Cuillandre, & Mellier 



Heckman 1996 ). Finally, even if our scenario turns 
out not to be applicable to spiral-dominated groups, 
it remains that ellipticals must eject gas in strong 
winds. Thus, our model may still be valid for clusters, 
which are elliptical-dominated. This will be explored 
in a subsequent work. 

We note as well that in our scenario, just as there is 
typically a large outflow from the halo, there is also a 
strong evaporative wind that ejects material from the 
Local Group. As a result, most baryons eventually 
reside in hot, intergalactic (as opposed to intragroup) 
gas. If this gas stays hot, it could perhaps be the 
ionized intergalactic (as suggested by Gunn-Peterson 
limits on the neutral intergalactic medium). If it does 
cool, it presents serious problems, as it would lead to 
prodigious but unobserved absorption of extragalactic 
radiation. 

Finally, we reiterate that stellar remnants and their 
associated hot ejecta are conservative candidates for 
both the halo microlensing objects and for the bary- 
onic dark matter. If these can be ruled out, then we 
are forced to conclude that the microlensing objects 
and the dark baryons are something stranger still. 
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FIGURE CAPTIONS 



1. Schematic diagram of model features. 

2. The total galactic star formation rate i/j = 
ncipcioud (solid curve) and number of clouds 
Uc (dashed curve) as a function of time. Two 
bursts of star formation are evident, one shortly 
after decoupling when the density is high, and 
one during halo collapse (5 Gyr). 

3. Mass evolution for (a) the halo, (b) the local 
Group. 

4. Metallicity evolution for the halo and Local 
Group. 

5. Halo metallicity distribution. The points are 



for globular clusters (Pagel 1988); only the low 
metallicity members {[Z] < —1) are included in 
the analysis. The theory (histogram) is binned 
in the same way as the data, and normalized to 
minimize x^. 

6. The halo white dwarf luminosity function, with 
L in units of Lq. Only results for halo white 
dwarfs are shown. The two peaks correspond 
to the two bursts of star formation; the more 
luminous peak is the more recent burst at t = 5 
Gyr. This second peak has My > 20, well below 
current observational limits. 



This 2-column preprint was prepared with the AAS lATJiiX 
macros v4.0. 
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